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ABSTRACT 

We use the integrated polarized radio emission at 1.4 GHz (Hj, 4 GHz) from a large sample of AGN 
(796 sources at redshifts z < 0.7) to study the large-scale magnetic field properties of radio galaxies 
in relation to the host galaxy accretion state. We find a fundamental difference in Hi. 4 ghz between 
radiative-nrode AGN (i.e. high-excitation radio galaxies, HERGs, and radio-loud QSOs) and jet-mode 
AGN (i.e. low-excitation radio galaxies, LERGs). While LERGs can achieve a wide range of Hi. 4 ghz 
( up to ~30%), the HERGs and radio-loud QSOs are limited to n 14 GHz 15%. A difference in 
IIi .4 ghz is also seen when the sample is divided at 0.5% of the total Eddington-scaled accretion rate, 
where the weakly accreting sources can attain higher values of Hi. 4 ghz- We do not find any clear 
evidence that this is driven by intrinsic magnetic field differences of the different radio morphological 
classes. Instead, we attribute the differences in Hi. 4 ghz to the local environments of the radio sources, 
in terms of both the ambient gas density and the magnetoionic properties of this gas. Thus, not only 
are different large-scale gaseous environments potentially responsible for the different accretion states 
of HERGs and LERGs, we argue that the large-scale magnetised environments may also be important 
for the formation of powerful AGN jets. Upcoming high angular resolution and broadband radio 
polarization surveys will provide the high precision Faraday rotation measure and depolarization data 
required to robustly test this claim. 

Subject headings: radio continuum: galaxies - galaxies: active - galaxies: magnetic fields 


1. INTRODUCTION 

The interaction of powerful relativistic jets in radio- 
loud active galactic nuclei (AGN) with their environment 
is observed to strongly affect the dynamics of interstellar 
and intergalactic gas (e.g. Birzan et al. 2008), likely im¬ 
pacting the formation and evolution of the most massive 
galaxies in the Universe (e.g. Best et al. 2006; Croton et 
al. 2006). There exists a wide variety of radio-loud AGN 
jet types and powers, whose appearance is considered to 
be directly related to some combination of the intrinsic 
supermassive black hole properties (i.e. mass and spin), 
the amount of material feeding the black hole, the avail¬ 
ability of large amounts of magnetic flux, and the prop¬ 
erties of the environment into which the jet propagates 
(e.g. Stawarz et al. 2008, Laing & Bridle 2014, Zamani- 
nasab et al. 2014). One of the key open questions in the 
study of AGN is how the production of powerful rela¬ 
tivistic jets in radio-loud AGN is related to the accretion 
system of the host galaxy. While substantial progress 
has been made in our understanding of the relation be- 
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tween the nuclear emission properties and the large-scale 
radio emission in AGN (e.g. Rawlings & Saunders 1991, 
Laing et al. 1994, Hardcastle et al. 2007, Buttiglione et 
al. 2010, de Gasperin et al. 2011, van Velzen et al. 2013, 
Mingo et al. 2014), these studies have not yet considered 
the magnetic properties of the source or its environment, 
even though the magnetic field is well accepted as a key 
parameter in AGN physics (e.g. Meier 2002, Sikora & 
Begelman 2013). 

The host galaxies of radio-loud AGN display a large 
range of optical spectral signatures that can be grouped 
into at least two main modes (Hine & Longair 1979, 
Laing et al. 1994): High Excitation Radio Galaxies 
(HERGs) and Low Excitation Radio Galaxies (LERGs). 
The HERGs have strong, high-ionisation emission lines 
and represent the classical picture of an AGN with an 
optically-thick, geometrically-thin, radiatively-efficient 
accretion disk (Shakura & Sunyaev 1973). This class 
represents the most optically-luminous AGN (including 
radio-loud QSOs), a population that can be easily iden¬ 
tified out to high redshift (z > 3). This class of AGN are 
also sometimes referred to as ‘radiative-mode’, ‘quasar¬ 
mode’ or ‘cold-mode’ AGN. On the other hand, in the 
LERG class of objects, the optical emission lines are 
very weak or non-existent, but the presence of radio jets 
clearly identify the source as an AGN. In this case, the 
presence of a radiatively-inefficient accretion flow is in¬ 
ferred (e.g. Narayan & Yi 1995). LERGs are also known 
as ‘jet-mode’, ‘radio-mode’ or ‘hot-mode’ AGN. HERGs 
and LERGs have been proposed to be fuelled in funda¬ 
mentally different ways (e.g. Best et al. 2005). The high 
radiative luminosities of HERGs are considered to be fu¬ 
elled by a large supply of cold gas for accretion, poten¬ 
tially provided by a merger with a gas-rich disk galaxy 
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or internal dynamical processes in the host galaxy it¬ 
self. The radiatively-inefficient, low-luminosity LERGs 
require much lower gas accretion rates, potentially pro¬ 
vided by Bondi accretion of the hot X-ray halos sur¬ 
rounding their massive elliptical galaxy hosts (Hardcastle 
et al. 2007). 

Jets from both HERGs and LERGs are generally con¬ 
sidered to be launched in similar manners from the vicin¬ 
ity of the central super-massive black hole of their host 
galaxy by large-scale, ordered magnetic fields (Bland- 
ford & Znajek 1977, Blandford & Payne 1982, Sikora & 
Begelman 2013). These intrinsically magnetised jets can 
propagate to large (^> kpc) distances, carrying energy 
and magnetic fields into their ambient galactic and ex- 
tragalactic environments. However, as the jets propagate 
through their host galaxy, the radio emission structure 
changes, presumably in part due to the type of inter¬ 
action with the ambient environment. In general, the 
large scale radio morphology is typically divided into two 
main types (Fanaroff & Riley 1974): low radio power 
Fanaroff-Riley Class 1 (FR1) sources where the inner jet 
and core have the brightest radio emission, and high ra¬ 
dio power FR2 sources where the brightest radio emis¬ 
sion is located at the jet termination ‘hotspots’. The 
traditional FR1 /FR2 radio luminosity density divide oc¬ 
curs at ~ 10 25 W Hz -1 (Fanaroff & Riley 1974), with 
a later study by Ledlow & Owen (1996) showing that 
the FR1/FR2 divide also depends on the host galaxy 
optical magnitude, with the division occurring at higher 
radio luminosities in more optically luminous hosts. This 
result indicates that the FR1/FR2 morphology is most 
likely due to a combination of jet power and environment, 
since it is thought that the Ledlow & Owen (1996) divi¬ 
sion arises due to jets being more easily disrupted into 
FR1 morphologies in the more massive host galaxies. 

The FR1 and FR2 populations also appear to differ in 
their accretion states, with HERGs often associated with 
FR2s and LERGs mainly of the FR1 type (e.g. Laing et 
al. 1994). The common interpretation of this association 
is that with a higher accretion rate, due to the plentiful 
supply of cold gas, the HERGs can produce more pow¬ 
erful jets, that are more likely to develop into FR2-type 
sources; while the weakly accreting LERGs produce the 
less powerful FR1 jets. However, at low redshift, where 
both LERGs and HERGs can be equally well identi¬ 
fied, there are several examples of FRl-HERGs and FR2- 
LERGs (Hardcastle et al. 2007, Best 2009, Buttiglione et 
al. 2010). A detailed comparison of the FR.-type and host 
galaxy accretion mode found that the FR class was inde¬ 
pendent of the HERG/LERG divide (Gendre et al. 2013). 
Thus, the nature of the relationship between the accre¬ 
tion mode of the host galaxy, the formation of powerful 
radio jets, the morphological appearance of these jets on 
large scales, and the large scale environment is still un¬ 
clear. 

The radio polarization properties of radio-loud AGN 
are important because they probe the internal structure 
of the jet magnetic field, as well as providing a unique 
diagnostic of the magnetoionic environment of the radio 
source through Faraday rotation of the linear polariza¬ 
tion (e.g. Laing 1988, Garrington et al. 1988, Laing et 
al. 2008). The relationship between the large scale po¬ 
larization properties of the radio jets/lobes and the host 


galaxy accretion state is important not only because it is 
a sensitive probe of the local environment, and thus the 
origin of the gas fuelling the black hole, but it also probes 
the magnetic field properties of the gas, and therefore 
helps provide an assessment of the ability of the black 
hole to accumulate the sufficiently large amounts of mag¬ 
netic flux required for launching powerful radio-loud jets 
(e.g. Tchekhovskoy et al. 2012). 

In this paper, we present an analysis of the accretion 
properties of host galaxies (as probed by nuclear opti¬ 
cal emission lines) in relation to the large-scale mag¬ 
netic field properties of their AGN jet (as probed by 
linearly-polarized radio emission). In Section 2, we de¬ 
scribe the construction of our sample of polarized HERG 
and LERG sources, the classification of their total in¬ 
tensity radio morphologies and the measurement of their 
linear sizes. In Section 3, we present our results on the in¬ 
tegrated 1.4 GHz polarization properties of HERGs and 
LERGs. We also investigate the dependence of the polar¬ 
ization on radio morphology, linear size, Faraday rotation 
and the Eddington-scaled accretion rate. In Section 4, 
we discuss the implications of our results on the link be¬ 
tween the large scale magnetised environment of HERGs 
and LERGs and their black hole accretion state. Our 
conclusions are listed in Section 5. Throughout this pa¬ 
per, we assume a flat ACDM cosmology with Ho = 67.3 
km s -1 Mpc -1 , 9 m = 0.315 and Ha = 0.685 (Planck 
Collaboration et al. 2014). 

2. SAMPLE CONSTRUCTION 

The only lrigh-resolution radio survey covering a sig¬ 
nificant fraction of the sky (~82%) with polarization 
sensitivity is the NRAO VLA Sky Survey (NVSS), con¬ 
ducted at 1.4 GHz (Condon et al. 1998). The NVSS has 
an angular resolution of 45” with an approximately uni¬ 
form rms sensitivity of ^0.45 mJy beam -1 in Stokes I 
and ~0.29 mJy beam -1 in Stokes Q and U. While the 
NVSS catalog already lists the basic integrated polariza¬ 
tion properties, we have generated our own catalog by 
downloading the full NVSS IQU mosaics 8 and using the 
source-finding program AEGEAN (Hancock et al. 2012). 
This was done as part of a larger project in order to bet¬ 
ter characterise the morphologies and multi-component 
nature of polarized sources in the NVSS and has gener¬ 
ated catalog positions for components in polarized inten¬ 
sity, Stokes Q and Stokes U, independent of Stokes I. 
AEGEAN is an optimised source-finder designed for reli¬ 
ably identifying compact radio sources based on a Lapla- 
cian kernel. It finds peaks in an image based on a ‘seed’ 
threshold (cr s ) and grows an island of pixels around that 
peak based on a ‘flood’ threshold (erf), where the local 
image noise level (cr rms ) is estimated over a 20 x 20 beam 
area. A curvature image is then used to determine how 
many Gaussian components are used to describe the is¬ 
land, returning a catalog of component flux densities and 
angular sizes. For our polarization source finding in the 
NVSS polarized intensity images (p = ( Q 2 + U 2 ) 1 / 2 ), we 
used cjf = 7tr rms and a s = 8<r rms , with (j rms ~ 0.2 mJy 
being typical for most fields. This has resulted in the in¬ 
clusion of only those polarized sources with peak polar¬ 
ized intensities above 8a of the local noise level in polar¬ 
ized intensity. Polarization bias was corrected for using 

8 http://www.cv.nrao.edu/nvss/ 
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the estimator p 0 = ( p 2 — where po is an esti¬ 

mate of the true polarized intensity (Simmons & Stewart 
1985). This has produced catalogs of ^87,000 compo¬ 
nents in polarized intensity and total intensity (Stokes 
I), as well as Stokes Q and U individually, over the en¬ 
tire NVSS area. A subset of that catalog is used here. 
The full polarization catalog will be presented elsewhere. 

2.1. HERG/LERG parent sample 

In order to calculate physically interesting quantities, 
we need to match each polarized radio source with its 
spectroscopically identified optical host galaxy. The 
largest optical spectroscopic survey overlapping with the 
NVSS is the Sloan Digital Sky Survey (SDSS), which 
covers ~10,000 square degrees of the Northern Hemi¬ 
sphere sky (DR 7; Abazajian et al. 2009). In addition 
to this, the Faint Images of the Radio Sky at Twenty 
centimetres (FIRST) survey (Becker, White & Helfand 
1995), which overlaps with the SDSS survey area, has a 
similar sensitivity limit as the NVSS but a significantly 
higher angular resolution of 5”, that aids substantially 
in obtaining reliable cross-matches with the SDSS. Com¬ 
bining the NVSS, FIRST and SDSS catalogs, Best & 
Heckman (2012), hereafter BH12, produced a large sam¬ 
ple of radio-loud AGN, with spectroscopic classification 
into HERGs and LERGs out to a redshift z ~ 0.7. They 
used the emission line information from the value-added 
spectroscopic catalogs available at http://www.mpa- 
garching.mpg.de/SDSS/ (cf. Brinchmann et al. 2004) to 
robustly classify 481 HERGs and 9,863 LERGs out of 
18,286 radio-loud AGN. We also use the line emission 
information from these catalogs in this paper. 

2.2. HERG/LERG polarization sample 

We constructed the sample used in this paper by cross¬ 
matching our NVSS polarization catalog with the BH12 
catalog of HERGs and LERGs. Since BH12 have al¬ 
ready carefully cross-matched the NVSS sources with 
their SDSS optical host galaxies using the additional, 
higher-resolution information provided by FIRST, our 
catalog is immediately highly reliable. Our initial cat¬ 
alog used a cross-matching radius of 1 arcminute and 
produced 802 matches. After manual visual inspection 
of the optical (SDSS) and radio (NVSS & FIRST) im¬ 
ages of all matches, we found only 6 matches that had ei¬ 
ther misidentified optical hosts or were in confused fields 
(i.e. there was a potential contribution to Hi. 4 ghz by ra¬ 
dio emission from an unrelated radio source). Thus, the 
final catalog has 796 polarized radio sources, the vast 
majority of which are unresolved in the NVSS (~96%). 

In the case of multiple NVSS components related to 
one radio source, we first integrated the contribution of 
Stokes Q and U from each component before calcula¬ 
tion of the polarized intensity, and then divided by the 
sum of the total intensity components to get the inte¬ 
grated degree of polarization, Hi. 4 ghz- The polarization 
bias correction was applied to the integrated polarized 
intensity. Ideally the bias correction should be applied 
to each component first, but that would have prevented 
us from comparing the polarization properties of the re¬ 
solved sources with the unresolved sources in our analy¬ 
sis. In total, there are 32 resolved sources (in Stokes I) in 
our sample. Only 3 of these resolved sources have multi¬ 
ple polarized components, 19 have polarized component 


major axes greater than 45”, while the remainder have 
compact polarized components (major axis < 45”). For 
the 3 multiple-component polarized sources, the higher 
noise level of ctqu for the bias correction is calculated 
as the quadrature sum of the <tq and ajj noise for each 
component. 

Figure [l] shows examples of two SDSS galaxies (one 
HERG ana one LERG, with their accompanying optical 
spectrum shown) overlaid with the corresponding NVSS 
I and p images, as well as the FIRST total intensity im¬ 
age. Note how the radio images appear almost identical 
but the optical spectra differ considerably. 

2.3. FIRST morphological analysis 

The high angular resolution FIRST images allow us to 
classify the total intensity morphology of all the polar¬ 
ized sources, and to calculate the projected linear size 
from the total angular extent of the FIRST total inten¬ 
sity emission. We manually classified the radio morphol¬ 
ogy of the 796 polarized sources into FR1 and FR2 based 
on whether the source was edge-brightened with distinct 
hotspots (FR2) or edge-darkened with the highest sur¬ 
face brightness features along the jet and/or the core 
(FR1). To ensure that our results are not influenced 
by misclassifications, we denote source morphologies as 
‘FR1?’ or ‘FR2?’ in the cases where the FIRST image 
fidelity is not of sufficient quality to conclusively decide 
whether it really is an FR1 or FR2. We assign an FRO 
classification (e.g. Sadler et al. 2014) to those sources 
that are unresolved at 5”. This corresponds to a pro¬ 
jected linear extent of ~15 kpc at the median redshift 
of our sample (^median = 0.2). Thus, these sources are 
most likely compact symmetric objects (CSOs), blazars 
or sources with weak jets (e.g. low-luminosity AGN). The 
sources classified as ‘FRO?’ are ones that have a bright 
core but for which it is unclear if there is extended emis¬ 
sion that has been resolved out or is too faint to be 
detected in FIRST. In addition to the FR0/FR1/FR2 
classification, we also manually classified the radio mor¬ 
phology into ‘straight’, ‘bent’ and ‘compact’. Sources 
were determined to be ‘straight’ or ‘bent’ depending on 
whether the extended emission on either side of the host 
galaxy was co-linear or not. The criteria for the ‘com¬ 
pact’ classification is identical to the FRO classification. 
To calculate the projected linear size of the source in 
kpc, we use a manual estimate of the angular size of each 
source along with the redshift from the SDSS. When the 
angular size of a source was larger than the synthesised 
beam of the FIRST image, we estimated the full angular 
extent of the emission by drawing a straight line between 
the extremities of the emission. Unresolved sources were 
assigned an angular size of 5”, which corresponds to an 
upper limit on the linear size of the source. 

2.4. Infrared (IR) data 

We also cross-matched our sample of polarized HERGs 
and LERGs with the Wide-field Infrared Survey Explorer 
(WISE; Wright et al. 2010), conducted at four bands, 

3.4, 4.6, 12 and 22 /im, with 5<r point source sensitivities 
in unconfused regions of ^0.08, 0.11, 1.0, and 6.0 nrJy, 
and angular resolutions of 6.1, 6.4, 6.5 and 12 arcsec, 
in each band respectively. The IR properties across the 
WISE band are useful in determining the host galaxy 
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Fig. 1.— Example images from the cross-matched datasets, 
with the corresponding SDSS spectrum of the host galaxy 
included (obtained from the SDSS Science Archive Server, 
http://drl0.sdss3.org/). Top: Low Excitation Radio Galaxy 
(LERG) at (RA, Dec) = (172.52602, 19.74556). Bottom: High 
Excitation Radio Galaxy (HERG) at (RA, Dec) = (255.78027, 
24.62987). SDSS i-band (greyscale), NVSS total intensity (red 
contours), NVSS polarized intensity (green contours), FIRST total 
intensity (blue contours). 


types of our sample (e.g. quiescent, star-forming). We 
found 456 reliable cross-matches out of 796 sources down 
to a 5cr detection level in at least one of the four WISE 
bands. For conversion from the calibrated WISE VEGA 
magnitudes to flux density (in Jy), we followed the in¬ 
structions on the WISE webpage , 9 including the 10% flux 
correction in the 22 fj ,m band. 

3. RESULTS 

The integrated degree of polarization measured at 
1.4 GHz (Hi 4 GHz) in the NVSS can be influenced by 
a number of effects. Firstly, bright and resolved struc¬ 
tures in radio galaxies (i.e. knots, filaments, hotspots) 
typically do not achieve degrees of polarization larger 
than ~30% to 50% (Saikia & Salter 1988). This im¬ 
plies at least some amount of intrinsic disorder to the 
magnetic field on small scales in order to reduce it from 
the theoretical maximum of ~72% expected for optically 
thin synchrotron radiation in a completely uniform mag¬ 
netic field (e.g. Pacholczyk 1970). 10 Secondly, even with 
high intrinsic degrees of polarization on small scales, if 
the magnetic field structure is not globally uniform then 
polarization angle cancellation across the jet and lobe 
structure will significantly reduce the integrated degree 
of polarization. Thirdly, polarization measurements at 
1.4 GHz are expected to be strongly affected by Fara¬ 
day rotation, since the rotation of the polarization angle 
depends on wavelength-squared. Large spatial fluctua¬ 
tions in Faraday rotation are observed in radio galaxies 
(e.g. Laing et al. 2008), that cause significant depolariza¬ 
tion at 1.4 GHz. Finally, even if the intrinsic magnetic 
field is globally uniform and the Faraday rotation is negli¬ 
gible, observed asymmetries of the jet and lobe emission 
structure, due to physical bends or relativistic effects, 
will lead to polarization angle cancellation that can re¬ 
duce the degree of polarization when integrated over the 
entire source. 

For the 796 polarized sources, we find that Hj 4 GHz 
ranges from > 0.1% to < 30% (Fig. [2]), with a median 
value of 6.2%. This is consistent with the maximum 
integrated polarization expected at 1.4 GHz of < 30% 
from surveys of extragalactic radio sources (e.g. Hales 
et al. 2014, although a small number of exceptions were 
presented by Shi et al. 2010). Since our only available 
polarization measurements are at a single frequency and 
they have an angular resolution of 45” (corresponding 
to ~150 kpc at the median redshift of 0.2 of the sam¬ 
ple), we cannot study the depolarization characteristics 
of our sample or the polarization morphology in any de¬ 
tail. Thus, we attempt to gain insight into the nature 
of HERG and LERG sources through their integrated 
degree of polarization properties at 1.4 GHz, in conjunc¬ 
tion with their optical spectroscopic properties and total 
intensity radio morphology. 

3.1. Integrated Degree of Polarization of HERGs and 

LERGs 

Out of the 9,863 LERGs and 481 HERGs identified 
by BH12 out to z ~ 0.7, we detect 741 LERGs and 55 
HERGs with integrated linear polarized intensity at 1.4 

9 http://wise2.ipac.caltech.edu/docs/release/allsky/expsup/ 

10 Note that a completely disordered field that is compressed 
can also produce a high degree of polarization (e.g. Laing 1981). 
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GHz greater than 8 a (~1.6 niJy). There are 67 LERGs 
and 2 HERGs that have an integrated polarized flux den¬ 
sity less than 1.6 mJy. These are still 8 a detections but 
have lower local noise estimates than are typical for the 
entire NVSS. The faintest total intensity source we detect 
in polarization has a Stokes I flux density of 8.6 mJy. At 
this flux level, there are 393 HERGs and 7,661 LERGs 
meaning that ~90% of sources down to this flux level 
remain undetected in polarization. The integrated to¬ 
tal and polarized flux densities of the detected polarized 
sources are shown in Figure [2] with HERGs identified as 
blue square symbols and LERGs as red plus symbols; 
the dashed lines are used to highlight that the majority 
of polarized sources have integrated 1.4 GHz percentage 
polarizations (Hi. 4 ghz) between 1 and 30%. Both types 
of sources cover similar ranges in polarized intensity and 
the high fractional polarization sources are not all clus¬ 
tered at the detection limit (Fig. [2]). The rms residual 
instrumental polarization of the NvSS is estimated as 
~0.3% (Condon et al. 1998, Stil et al. 2014). Two LERGs 
and two HERGs have Hi. 4 ghz < 0.3% (Fig. 2), and 
we exclude them from our subsequent analysis. Thus, 
the remaining sample is composed of 53 HERGs and 739 
LERGs. 

The LERG sources span the full range of expected 
ni. 4 GHz from >0.3% to ~30%, while the detected 
HERGs are restricted to Hi. 4 ghz < 15% (Fig. [3]). The 
median Hi. 4 ghz for LERGs is 6.2% and 4.2% for HERGs. 
A two-sided Kolmogorov-Smirnov (KS) test of Hi. 4 ghz 
for HERGs and LERGs gives a probability or ‘p-value’ of 
0.4%, indicating a significant rejection, at approximately 
2.9cr, of the null hypothesis (i.e. the two datasets are not 
likely drawn from the same underlying distribution ). 11 
In the local universe, LERGs dominate over HERGs at 
low radio luminosity (BH12), meaning that LERGs are 
more numerous than HERGs at low radio flux density in 
our volume limited sample (Fig. [2]). Ricean bias effects 
the measurement of polarized intensity more strongly at 
low flux density (see Section 2). Therefore, one might 
expect the LERGs to be more strongly affected by the 
polarization bias than the HERGs. If we restrict our 
sample to those sources with polarized intensity greater 
than 8 (tq[/, then we effectively eliminate the effect of the 
polarization bias (e.g. Simmons & Stewart 1985). We 
also remove the 32 extended sources since the bias cor¬ 
rection applied to them is not ideal for resolved sources. 
This reduces the sample to 43 HERGs and 474 LERGs, 
with a KS test of Hi. 4 ghz for these HERGs and LERGs 
now giving a p-value of 0.2% (^3.la). This small in¬ 
crease in significance shows that the polarization bias 
does not strongly affect our full sample. 

A more important bias is the selection bias introduced 
by the threshold in polarized intensity at ~1.6 mJy (see 
Fig. 2). This causes the median Hi 4 GHz to increase for 
fainter sources because only sources with high values of 
ni. 4 GHz will be detected. If we consider the sources 
for which this threshold in polarized intensity is not 

11 Throughout this paper we use the two-sample KS test to de¬ 
termine the probability (or ‘p-value’) that the null hypothesis can 
be rejected for two samples being tested. In all cases, the null 
hypothesis is that the two samples are drawn from the same un¬ 
derlying distribution. In addition to the p-value, we also quote the 
significance level in terms of the equivalent result from a normally 
distributed process. 


as important (e.g. for sources brighter than 100 mJy), 
then we find no statistically significant difference be¬ 
tween Hi 4 GHz for HERGs and LERGs (see Appendix 
A). More sensitive polarization observations are required 
(by at least an order of magnitude) to conclusively test 
for a significant difference in n x 4 GHz between HERGs 
and LERGs down to I ~ 10 mJy. However, Stil et 
al. (2014) circumvented this problem somewhat by stack¬ 
ing NVSS sources in polarized intensity to obtain the 
median degree of polarization for sources that were too 
faint to be detected. They found that the median de¬ 
gree of polarization increases with decreasing flux den¬ 
sity (ni. 4 GHz OC 5^)2 ghz, where 5 i. 4 ghz is the NVSS 
flux density at 1.4 GHz). Our results suggest that one of 
the most likely causes of this increase in Hi.4ghz is the 
large number of LERGs at low flux with Hi. 4 ghz > 15%. 
Thus, the main focus of this paper is on the absence of 
HERGs with Hi. 4 ghz > 15% and the ability of some 
LERGs to achieve values of Hi.4ghz up to ~30%. 

3.1.1. Integrated degree of polarization of radio-loud QSOs 

Hammond et al. (2012) presented a catalog of 4,003 
polarized radio sources with redshifts, 815 of which were 
identified as QSOs in the SDSS, with 89 at 2 < 0.5. 
Radio-loud QSOs are the equivalent of HERGs (ie. radio- 
loud AGN with a radiatively-efficient accretion disk) but 
with a more direct line of sight to the AGN nucleus 
and typically more luminous in the optical, allowing 
them to be more easily detected out to high redshift 
(e.g. van Velzen et al. 2015). A clear difference be¬ 
tween the 1.4 GHz fractional polarization distributions 
of SDSS galaxies and of QSOs was found by Hammond 
et al. (2012), with the SDSS QSOs not exceeding frac¬ 
tional polarizations greater than ~ 15% while the SDSS 
galaxies could reach as high as ~ 30%. Here we have 
essentially expanded on their sample of SDSS galaxies 
and separated them into the more physically meaning¬ 
ful classes of HERGs and LERGs. In Figure [3j sepa¬ 
rate histograms of Hi.4ghz are shown for the LERGs, 
HERGs and QSOs. Now there are roughly equal num¬ 
bers of radiative-mode AGN (i.e. HERGs and QSOs) and 
jet-mode AGN (i.e. LERGs), and the maximum degree of 
polarization of the radiative-mode AGN is still limited to 
< 15%. This strongly suggests a fundamental difference 
between the polarization properties of radiative-mode 
AGN and jet-mode AGN at 1.4 GHz. We note that there 
are potentially strong redshift and evolutionary effects 
for the radio-loud QSO sample (0.06 < 2 < 5.3). Thus, 
we limit most of our detailed analysis and discussion to 
the polarized HERGs and LERGs because they are sim¬ 
ilarly distributed in redshift, with 2 me dian,HERG = 0.25 
and 2 me dian,LERG = 0.18 (Figured). In Appendix A, we 
replot Figure [I] for the sources whh Stokes I > 100 mJy 
only. 

3.2. Radio Luminosity Distribution of Polarized 
HERGs and LERGs 

The distribution of n^GHz with 1.4 GHz radio lu¬ 
minosity (L i.4 ghz) is shown in Figure [5] We calculate 
the radio luminosity using Li. 4 ghz = 47 tD| j S'i. 4 ghz (1 + 
2 ) _ (“ +1 ), where Dl is the luminosity distance, Si. 4 GHz 
is the integrated total intensity and we have assumed 
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Fig. 2.— Integrated linear polarized flux density ( p ) versus to¬ 
tal flux density (/) for HERGs (blue squares) and LERGs (red 
plus symbols). Diagonal dashed lines represent constant values of 
integrated degree of polarization. 
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Fig. 3.— Histograms of the integrated degree of polarization 
(Ili. 4 GHz) for LERGs (red), HERGs (blue, hatched) and radio- 
loud QSOs (grey). 
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Fig. 4. — Integrated degree of polarization at 1.4 GHz (ni 4GHz) 
versus redshift (z) for LERGs (red), HERGs (blue, hatched) and 
radio-loud QSOs (grey). 

a spectral index a = —0.7. 12 In the local universe, 
LERGs dominate the radio-loud AGN number counts at 
low luminosities (22 < log(Li. 4 GHz [W Hz -1 ]) < 25), 
while above Li. 4 GHz ~ 10 26 W Hz -1 , HERGs begin 
to dominate (BH12). While examples of both popu¬ 
lations are found across the full range of HERG and 
LERG luminosities, we only detect polarized HERGs 
above Li. 4 GHz ~ 10 24 W Hz -1 . 

Since LERGs are more numerous than HERGs at lower 
luminosity and we have just shown that LERGs can 

12 We define the spectral index, a, such that the observed flux 
density (S) at frequency v follows the relation S v oc iA“. 


Fig. 5.- Integrated degree of polarization at 1.4 GHz (Hi 4qhz) 
versus the 1.4 GHz radio luminosity, in units of W Hz -1 , for 
HERGs (blue squares) and LERGs (red plus symbols). Note that 
we are insensitive to many low-luminosity, 1 ow-Hi. 4 ghz sources 
due to the polarization detection limit. 


achieve higher values of Hi 4 ghz> then one might expect 
an anti-correlation between Hi. 4 ghz and Li. 4 GHz- We 
are insensitive to low-luminosity low-^ ,4 GHz sources due 
to the polarization detection limit, and this strongly af¬ 
fects any robust inferences on any dependence of Hi. 4 ghz 
on Li. 4 ghz- If we only consider sources with total inten¬ 
sities greater than 50 mJy, we have roughly 70% of both 
HERG and LERG sources below the detection limit in 
polarized flux density. This is roughly the limit at which 
useful statistics can be derived from data with upper 
limits (e.g. Antweiler & Taylor 2008). Using a Kendall’s 
tau rank correlation, we find a statistically significant 
(p-value of 6 x 10 -8 , ~5.4a), but weak, anti-correlation 
(correlation coefficient of — 0 . 1 ) between Hi. 4 ghz and 
Li. 4 ghz for all sources. Considering the HERG sources 
only, we find no statistically significant correlation but 
we find a marginally significant anti-correlation between 
ni. 4 GHz and Li. 4 GHz for the LERG sources (p-value of 
0.013%, ^3. 8 a, for a correlation coefficient of —0.1). 
There is no evidence for an anti-correlation between 
Hi .4 GHz and Li. 4 ghz for the brightest sources (i.e. with 
Stokes I > 100 mJy, as shown in Appendix A). This 
suggests that any anti-correlation between Hi. 4 ghz and 
Li. 4 GHz is mainly driven by the change in the domi¬ 
nant population of radio-loud AGN towards lower lu¬ 
minosities from HERGs to LERGs. Much deeper radio 
polarization observations are required to more robustly 
test the existence of an anti-correlation between Hj 4 GHz 
and Li. 4 ghz for the LERG sources. Of course, radio 
luminosity may be anti-correlated with Hi. 4 ghz inde¬ 
pendent of the HERG/LERG classification and deeper 
polarization observations would also provide a key test 
for the continued absence of low-luminosity HERGs with 
ni.4GHz > 15%. 

It is interesting to extend this plot to higher radio 
luminosities by including all the known radio polarized 
QSOs (Hammond et al. 2012) that overlap with the sky 
coverage of our sample (i.e. the SDSS area). Figure | 6 | 
shows Hi 4 GHz versus Li. 4 GHz again but with the QSOs 
included (grey squares). This demonstrates the exten¬ 
sion of the radiative-mode AGN population to higher lu¬ 
minosities and clearly shows the limitation of radiative¬ 
mode AGN to values of Hi. 4 ghz < 15%, across five or¬ 
ders of magnitude in radio luminosity. 
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Fig. 6.— Integrated degree of polarization at 1.4 GHz (ni. 4 GHz) 
versus the 1.4 GHz radio luminosity, in units of W Hz -1 , for 
HERGs (blue squares) and LERGs (red plus symbols) and radio- 
loud QSOs (grey squares). 




Fig. 7. — Integrated degree of polarization at 1.4 GHz (ni4GHz) 
versus the 1.4 GHz radio luminosity, in units of W Hz -1 , for FRO 
(grey triangles), FR1 (red circles) and FR2 (blue diamonds) radio 
morphologies. Filled symbols represent robust classifications, with 
open symbols representing uncertain classifications. See Sections 
2.3 and 3.3.1 for details. 



Fig. 8. — Empirical cumulative distribution functions (ECDFs) 
of the integrated degree of polarization (Hi. 4 ghz) 3 for FR2 (blue), 
FR1 (red) and FRO (grey) classifications. 

3.3. Radio morphology and integrated degree of 
polarization 

We now present an investigation of the origin of the 
difference in Hiaghz for HERGs and LERGs in relation 
to the morphological properties of the large scale radio 
emission (total intensity only, from the FIRST survey). 

3.3.1. FRO, FR1 and FR2 

Using our manual classification of FR0/FR1/FR2, de¬ 
scribed in Section 2.3, we find 14/6/18 HERGs that 



Fig. 9. — Empirical cumulative distribution functions (ECDFs) 
of the integrated degree of polarization at 1.4 GHz (Hi. 4 ghz) for 
radio morphologies classified as straight (green), bent (magenta) 
and compact (black). See Sections 2.3 and 3.3.2 for details. 



Fig. 10 . — Integrated degree of polarization at 1.4 GHz (Hi 4 GHz) 
versus projected linear source size in kpc, for LERGs (red plus 
symbols, with upper limits denoted by red triangles) and HERGs 
(blue squares, with upper limits denoted by blue triangles). See 
Sections 2.3 and 3.3.3 for details. 



Fig. 11.— Empirical cumulative distribution function (ECDF) of 
the source linear size in kpc (excluding sources with upper limits), 
for LERGs with Hi 4ghz < 15% (red solid line), HERGs (blue 
solid line) and LERGs with ni. 4 GHz > 15% (red dashed line). 


are FR0/FR1/FR2 with 15 uncertain classifications and 
89/221/174 LERGs that are FR0/FR1/FR2 with 255 
having an uncertain classification. It is immediately clear 
from this that the HERG/LERG class is essentially in¬ 
dependent of the FR type. This is consistent with other 
studies by Best et al. (2009), Lin et al. (2010) and Gen- 
dre et al. (2013). Figure [7] again plots Hi. 4 ghz versus 
radio luminosity but now with the FR-type identified. It 
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is clear from this plot that the range of IIi.4ghz is not 
restricted to values less than 15% for FR1, FR2 or FRO 
sources. Thus, the HERG-LERG difference in IIi.4ghz 
cannot be easily explained by the radio-morphology type. 
Similar results are obtained if we only consider the 
brightest sources, with I > 100 mJy (see Appendix 
A). Figure [ 8 ] shows the empirical cumulative distribu¬ 
tion functions (ECDFs) for IIi.4GHz of FRO, FR1 and 
FR2 sources. Considering only the definitive FR1 and 
FR2 type sources, we find a marginally significant dif¬ 
ference (~2.3<r) in IR 4 GHz between the two types, with 
FRls having a median IR .4ghz,fri of 6.4% and a median 
IIi. 4 GHz,FR 2 of 5.3% for FR2s. Similar results were pre¬ 
viously found by Grant (2011) for a smaller number of 
sources (19 FRls and 17 FR2s). We find no significant 
difference in IR.4ghz between FRO and FR1 sources (p- 
value of 18%) but get a ~3<r difference (p-value of 0.3%) 
in IR 4 GHz between FRO and FR2 sources. The median 
IIi. 4GHz ,fro of FRO sources is 6 . 1 %. 

3.3.2. Straight jets and bent jets 

Intrinsic polarization angle cancellation in radio 
sources with bent or disturbed morphologies could result 
in low integrated degrees of polarization, while the undis¬ 
turbed, straight jet sources may achieve much higher val¬ 
ues of IIi. 4 GHz- Figure[9]presents ECDFs of IIi.4ghz split 
into straight, bent and compact sources (as defined in 
Section 2.3). The majority of polarized sources are clas¬ 
sified as straight (58%), with 25% bent and the remain¬ 
ing 17% being compact. For sources with IIi .4 ghz > 15% 
the fraction of straight sources is 60% (similar to the frac¬ 
tion for all sources), however, the fraction of bent sources 
falls to only 5%. KS-tests show that both compact and 
straight sources have significantly different IR. 4 ghz dis¬ 
tributions than bent sources (~4.7er and 4.3er, respec¬ 
tively) while compact and straight sources do not differ 
significantly from each other (p-value of 3%). We find 
that 40% of HERGs are straight, 32% are bent and 28% 
are compact, while 60% of LERGs are straight with 24% 
bent and 16% compact. 

3.3.3. Influence of the ‘radio core’ 

The contribution of radio emission from active or 
‘restarted’ emission regions near the central engine 
(i.e. the radio core) could decrease the integrated degree 
of polarization. This is possible because the inner re¬ 
gions of the radio source should suffer larger amounts of 
depolarization due to being more deeply embedded in the 
host galaxy. Furthermore, synchrotron self-absorbed and 
optically thick emission regions can only reach intrinsic 
degrees of polarization of approximately 10% (e.g. Pa- 
cholczyk 1970). In this case the additional total flux 
density provided by the core is likely to contribute very 
little additional polarized flux, causing the overall degree 
of polarization to decrease. To investigate whether or not 
this was an important effect for our sample, we deter¬ 
mined the number of extended sources (i.e. with angular 
sizes >10”) that also had radio cores in the FIRST im¬ 
ages. Out of the 38 polarized HERGs with angular sizes 
greater than 10 ”, 6 of them have radio cores (~16%). 
The median value of Hi. 4 ghz for these 6 sources is 3.7%, 
while the median value of Hi. 4 ghz for the other 32 ex¬ 
tended sources without a core is 4.4%. This indicates 


that the presence of a core may have a small effect in re¬ 
ducing Hi 4 ghz for HERGs (although a two-sided KS test 
finds no significant difference in Hi. 4 ghz for HERGs with 
and without a core, p-value of 58%). Of the 597 extended 
LERGs, 137 have radio cores (~23%). This shows that 
the presence of a core is slightly more likely for LERGs 
that HERGs. The median value of Hi. 4 ghz for LERGs 
with cores is 6 . 6 % and without cores is slightly lower, 
at 6.1%. In fact, there are 17 LERGs with radio cores 
that have Hi. 4 ghz > 15%, showing that the presence of 
a core is not a hindrance to high integrated degrees of 
polarization; these cores may be dominated by bright, 
highly polarized, optically-thin inner jet regions. 


3.3.4. Projected linear size 

One might expect the extended polarized emission of 
sources with larger linear sizes to have less Faraday depo¬ 
larization, and thus higher Hi. 4 ghz, due to the propaga¬ 
tion of the radiation through a presumably less dense 
magnetoionic environment (e.g. Hardcastle & Krause 
2014). To investigate this possibility, we compare the 
projected linear extent of the radio emission in the source 
rest frame, measured as described in Section 2.3, to 
III. 4 GHz, for the HERG and LERG sources (Figure [l0|. 
For all sources, there is no obvious strong dependence of 
III. 4 GHz on linear size. Excluding the sources with up¬ 
per limits, the median linear size of the HERG sample is 
102.8 kpc, and 82.8 kpc for the LERG sample. Thus, if 
the Hardcastle & Krause (2014) simulation results were 
true for all sources, then one would expect Hi. 4 ghz for 
HERGs to extend to higher values than the LERGs, com¬ 
pletely the opposite to what is observed. A two-sided KS 
test indicates, however, that the linear size distribution 
of the HERGs is not significantly different from that of 
the LERGs (p-value of 17%). The median linear size of 
sources with Hi. 4 ghz > 15% is 64.7 kpc, with a KS-test 
indicating a significant difference in linear size compared 
to sources with Hi. 4 ghz < 15% (p-value of 0.6%, ~2.7 a). 
This surprisingly indicates that the most highly polarized 
LERGs are the ones with smaller projected linear extents 
(Figure 11). If we include the upper limits as measure¬ 


ments then the median linear size for HERGs decreases 
to 75.4 kpc and 72.0 kpc for LERGs. The larger decrease 
for HERGs occurs because a greater fraction of HERG 
sources are compact (c.f. Section 3.3.2). 

We caution that the smaller linear size estimates in the 
fainter NVSS sources may be confounded by a smaller es¬ 
timate of the angular size due to the extended emission 
being resolved out/undetected in FIRST. Furthermore, 
while the measurement of Hi. 4 ghz is independent of the 
FIRST angular size measurement, it is worth consider¬ 
ing here whether or not the NVSS angular size influences 
the measurement of Hj 4 GHz differently for HERGs or 
LERGs, even though the majority of sources are unre¬ 
solved (84%). However, a KS test of NVSS angular size 
for HERGs and LERGs shows no significant difference 
(p-value of 26%). 


3.4. Environmental diagnostics of HERGs and LERGs 

Differences in the large scale gaseous and magnetoionic 
environments are prime candidates for explaining the dif¬ 
ference in Hi 4 ghz for HERGs and LERGs, since the 
radio jet emission and morphology is expected to be 
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strongly influenced by its local environment (e.g. Laing 
et al. 2014). Here we explore several potential tracers of 
the environment local to the source. 

3.4.1. Faraday Rotation Measure 

The integrated fractional polarization at 1.4 GHz is ex¬ 
pected to be strongly affected by Faraday rotation, where 
magnetoionic material along the line of sight causes a ro¬ 
tation of the intrinsic linear polarization angle (T) with 
wavelength squared, AT = RMA 2 . The Faraday rota¬ 
tion measure (RM) is defined as 0.81 n e B\\dl rad m~ 2 , 
where n e is the electron number density in cm -3 , B\ is 
the line-of-sight magnetic field strength in /xG and L is 
the distance through the magnetoionic region in parsecs. 
For example, with B\\ ~ 1 /xG, n e ~ 10~ 4 cm -3 and 
L ~ 100 kpc gives an RM of ~±8 rad m -2 , which causes 
a polarization angle rotation of 20 degrees at 1.4 GHz. 
Fluctuations in RM of this order (and larger) across 
the jets and lobes of radio galaxies are not uncommon 
(e.g. Laing et al. 2008, Guidetti et al. 2012). Thus, we 
expect such RM fluctuations local to the radio source to 
strongly reduce the integrated fractional polarization at 
1.4 GHz. 

In order to investigate the relative importance of this 
for HERGs and LERGs, we use the RM catalog of Taylor 
et al. (2009), that produced 37,543 RMs by reprocessing 
the NVSS data into two channels centred on 1.365 and 
1.435 GHz. Using this catalog, we can compare the RM 
of the integrated polarized emission for our sample of 
HERGs and LERGs. In Figure [T2j we plot the empirical 
cumulative distribution functions (ECDFs) of the abso¬ 
lute value of the RMs for 295 LERGs and 34 HERGs 13 . 
A two-sided KS test finds no significant difference be¬ 
tween the two distributions (p-value of 95%). Further¬ 
more, we also find no significant difference between the 
|RM| of sources with Hi. 4 gh z > 15% and those with 
ni. 4 GHz < 15% (p-value of 99%). We searched for differ¬ 
ences in |RM| between FR types, between straight and 
bent sources and in relation to linear size, finding no 
significant difference or dependence on any of the sub¬ 
samples. This is somewhat unsurprising since we know 
that the RMs from Taylor et al. (2009) are dominated 
by the Galactic contribution (Schnitzeler et al. 2010, 
Stil et al. 2011). Attempts to model the Galactic RM 
component and subtract it from the total RM to ro¬ 
bustly recover the extragalactic contribution is not possi¬ 
ble with the current limited data due to poor constraints 
on individual source errors (e.g. Oppermann et al. 2015). 
However, this is expected to be possible in future with 
higher precision RM measurements from observations 
with much broader and more continuous wavelength- 
squared coverage. 

3.4.2. Depolarization due to the local environment 

As described above, fluctuations in the magnetoionic 
medium across the source can significantly depolarise the 
radio emission. The recent study of Fames et al. (2014) 
presented a parametrization of the change in the degree 

13 The reason we did not find RMs for all our sources is that the 
Taylor et al. (2009) RM catalog uses a higher threshold of 8 oqij 
(~2.3 mjy). 



Fig. 12. — Empirical cumulative distribution function (ECDF) 
of the absolute value of the Faraday rotation measure (|RM|) in 
rad/m -2 , for LERGs (red solid line), HERGs (blue solid line) and 
LERGs with ni 4 GHz > 15% (red dashed line). 

of polarization with wavelength (i.e. a degree of polar¬ 
ization spectral index, /3, where oc X 13 ) for 951 ra¬ 
dio sources. They found that sources with a steep total 
intensity spectrum exhibited depolarization (a decreas¬ 
ing degree of polarization with increasing wavelength), 
while flat-spectrum sources typically maintain approxi¬ 
mately constant degrees of polarization over large ranges 
in wavelength. From this they concluded that /3 is pre¬ 
dominantly affected by the magnetoionic environment lo¬ 
cal to the source, rather than by material in the distant 
foreground (i.e. intergalactic magnetic fields and/or the 
Milky Way). Only twelve of the sources presented here 
have an estimate of j3 from Fames et al. (2014), with (3 
ranging from —0.8 to 0.3, with a median value of —0.2. 
While this is insufficient for a robust statistical analysis 
of our sample, we can use the overall results of Fames et 
al. (2014) to infer that Hi. 4 ghz for the majority of our 
sources, which have steep total intensity spectral indices, 
are likely affected by depolarization due to the local envi¬ 
ronment. Further observations are required to determine 
if HERGs and LERGs have different depolarisation prop¬ 
erties. 

3.4.3. Galaxy number density 

We now consider the polarization properties of HERGs 
and LERGs in comparison with one of the most direct 
measure of environment, galaxy number density. Esti¬ 
mates of the galaxy number density are available for 
~20% of our polarized sources, from Lin et al. (2010). 
They estimated the excess number of galaxies over the 
mean background within 0.5 Mpc (So.sMpc) of the host 
galaxies of 1040 extended FIRST sources. For the polar¬ 
ized HERGs and LERGs, the median value of E 0 . 5 m P c 
is 6.6, compared to the unpolarized sources, that have 
a median So .5 Mpc of 8.8, with a KS test indicating a 
significant difference (p-value of 4.8 x 10 -5 , ~4 a). This 
means that the polarized sources are typically in more 
underdense environments, in comparison to the unpo¬ 
larised sources. The median value of Eo. 5 Mpc for all (po¬ 
larized and unpolarized) HERGs is 4.1 compared to 8.6 
for all LERGs, supporting previous studies that found 
that HERGs are generally in less dense environments 
that LERGs (e.g. Best 2004). A KS test of E 0 . 5 m P c for 
polarized HERGs versus polarized LERGs find no sig¬ 
nificant difference between the sources (p-value of 10%). 
However, only five polarized HERGs have Eg .5 Mpc esti- 
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Fig. 13.— Integrated degree of polarization at 1.4 GHz 
(Hi. 4 GHz) versus the total intensity spectral index (a) taken from 
Fames et al. (2014), for HERGs (blue squares) and LERGs (red 
plus symbols). See Section 4.1.2 for details. 

mated, with all having values of Ho. 5 Mpc < 6 . 8 . For the 
LERGs, the median values of IIi. 4 ghz are 6.3% and 7.5% 
for sources with So.sMpc greater than and less than the 
median value of Eo. 5 Mpc, respectively (p-value of 0.3%, 
~2.9 ct significance). This suggests that, at least for the 
LERGs, the highest integrated degrees of polarization 
appear to favour the lowest galaxy density environments. 

3.4.4. Spectral index as an environment density probe 

Samples of radio galaxies selected by spectral index 
have previously shown that the steepest spectral index 
sources typically reside in rich clusters of galaxies (Bald¬ 
win & Scott 1973; Slee, Siegman & Wilson 1983). This 
has been interpreted as a result of pressure-confinement 
of the radio lobes propagating in dense environments 
(e.g. Klamer et al. 2006). We have total intensity spec¬ 
tral index measurements (a) from Fames et al. (2014) 
for 32% of LERGs and 56% of HERGs, with 95% of 
LERGs and 100% of HERGs having a within the range 
— 1.5 < a < —0.3. This supports the expectation that 
the majority of our polarized sources are indeed lobed 
radio galaxies with the dominant polarized emission pre¬ 
sumably coming from the steep spectrum, extended re¬ 
gions of the source (c.f. Banfield et al. 2011, Hales et 
al. 2014). 

From the plot of Hi. 4 ghz versus a, shown in Figure [T3] 
one can immediately notice the absence of sources with 
large Hi. 4 ghz and steep a (i.e. a < —0.7). However, 
this appears to be only true for the LERG sources since 
the HERGs essentially uniformly cover their Hi. 4 GHz _ a 
parameter space. Using a Kendall’s tau rank correla¬ 
tion test, we find a significant correlation between a and 
Hi .4 GHz for the steep-spectrum LERGs, with a correla¬ 
tion coefficient of 0.2 at ~4.5<r significance. Employing 
the same correlation test for the HERGs only, finds no 
significant correlation. There are no known selection ef¬ 
fects in Fames et al. (2014) that could connect both the 
spectral index (derived from multiple surveys), and the 
polarisation fraction (derived from just NVSS). There¬ 
fore, if very steep values of a are indicative of dense en¬ 
vironments, then it is interesting to relate the absence 
of highly polarized LERG sources with a < —0.7 to the 
effect of the environment. 

3.5. Host galaxy and accretion rate properties of 
polarized HERGs and LERGs 
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In order to infer general properties of the HERG and 
LERG populations from the subset of sources with signif¬ 
icant integrated polarized emission at 1.4 GHz, we need 
to assess whether or not the polarized sources are gen¬ 
erally representative of their individual classes. Previous 
studies have found that in the local universe HERGs are 
generally found in disk galaxies with ongoing star forma¬ 
tion and with a high density of cold gas accreting at a rate 
greater than ~1% of the Eddington rate. In contrast, 
LERGs are mainly found in massive elliptical galaxies 
with more massive black holes fuelled by slowly cooling 
hot gas that accretes at less than 1% of Eddington. The 
accretion energy of LERGs is mainly channeled into their 
jets, that in turn limit the amount of hot gas cooling in 
their host galaxy and cluster environments in a feedback 
cycle. See Heckman & Best (2014) for a comprehensive 
review. In the remaining part of this section, we present 
the results of our analysis of the optical, infrared and ra¬ 
dio properties of the polarized HERG/LERG sample in 
order to determine their host galaxy and accretion rate 
properties. 


3.5.1. WISE Colour-Colour Host Galaxy Diagnostics 

The IR properties across the WISE band can be used 
to help separate early-type elliptical galaxies from galax¬ 
ies with a substantial amount of on-going star forma¬ 
tion in disky/spiral galaxies as well as identifying those 
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AGN with a bright accretion disk (i.e. a QSO). Figure 
shows the WISE colour-colour plot for our sample. The 
dotted horizontal and vertical lines are from Wright et 
al. ( 2010 ), who divide elliptical and spiral galaxies at 
a WISE [4.6]- [12] colour of +1.5 mag and find that the 
most powerful optical AGN lie above a [3.4]—[4.6] colour 
of +0.6 mag. From this we can see that the majority 
of LERGs are hosted by elliptical galaxies with only a 
small fraction classified as having spiral-type hosts. The 
LERGs classified as QSOs are potentially heavily ob¬ 
scured Type 2 AGN. The polarized HERGs are mainly 
in the WISE late-type category (72%) with most of the 
rest (22%) clustered within 0.5 mag of the dividing line 
between WISE early and late-type categories. This is 
broadly consistent with the BH12 detailed study of the 
optical host galaxy properties, where HERG host galax¬ 
ies were typically of lower stellar mass, lower black hole 
mass (i.e. less evolved) and had bluer colours than the 
LERGs. 


3.5.2. WISE Mid-IR Luminosity 

Gurkan et al. (2014) analysed the IR properties of a 
large sample (346) of radio-loud AGN including both 
nearby and high redshift sources (0.003 < z < 3.4). 
They found that their sample could be clearly divided 
into HERGs and LERGs in the mid-IR (22/ira) radio¬ 
luminosity plane, with a critical mid-IR luminosity of 
L 22 nm ~ 5 x 10 43 erg s _1 . If we divide our sample 
(z < 0.7) based on this mid-IR luminosity, we find 
that 75% of LERGs have + 22 ^™, below 5 x 10 43 erg s ~ 4 
and 78% of HERGs have mid-IR luminosities above this 
value. Figure [15] shows there is a large scatter of 1 to 2 
orders of magnitude about this dividing luminosity, how¬ 
ever, it may be a potentially useful diagnostic tool in the 
absence of optical spectroscopy. 

3.5.3. Eddington-Scaled Accretion Rate 
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Fig. 14. — WISE infrared colour-colour plot for LERGs (red plus 
symbols) and HERGs (blue squares). The dashed lines denote the 
regions where different galaxy populations are expected to be. See 
Section 3.5.1 for more details. 



Fig. 15.— Mid-IR luminosity (7.2 2 , , in erg/s) versus radio lu¬ 

minosity (Li. 4 GHzj in erg/s), for LERGs (red plus symbols) and 
HERGs (blue squares). The dashed line corresponds to the critical 
mid-IR. luminosity of 5x 10 43 erg/s defined by Gurkan et al. (2014). 
See Section 3.5.2 for details. 

Recent studies have investigated the idea that a fun¬ 
damental switch in the type of radio-loud AGN activity 
occurs at an accretion rate of ~1% of the Eddington 
accretion rate (e.g. Mingo et al. 2014). At Eddington 
ratios greater than 1%, the accretion disk structure is 
expected to represent the classical picture of an AGN in 
unified models with a geometrically-thin accretion disk 
(e.g. Urry & Padovani 1995), while at accretion rates 
much less than 1% of Eddington there is a geometrically- 
thick radiatively-inefhcient accretion flow (e.g. Narayan 
& Yi 1995). 

In order to estimate the Eddington-scaled total accre¬ 
tion rate for the polarized HERGs and LERGs, we follow 
the approach adopted by BH12. They combined the opti¬ 
cal bolometric luminosity, estimated from the OIII emis¬ 
sion line (Lboi,om); with the total mechanical luminosity 
of the radio jet (-L mec h), normalised to the Eddington lu¬ 
minosity (LEdd)i t° obtain the Eddington-scaled total 
accretion rate as (L mec h + Tboi,oin)/iEdd- 

To derive a mean bolometric correction for the SDSS 
[OIII] 5007 emission line luminosity, Heckman et al. (2004) 
used a multi-wavelength analysis of a large sample of 
powerful Type 1 AGN to find 

Tboi.om ~ 3500 Loiii, ( 1 ) 

with a scatter in Tboi.oiii of ±0.38 dex. The [OIII] 5 oo 7 


emission line is used because it is the strongest emission 
line in the majority of SDSS AGNs, and has minimal 
contamination from star-formation. Since this relation 
was derived mainly at low redshift [z ~ 0.1) its applica¬ 
bility up to redshifts of z ~ 0.5 is questionable. However, 
Bernardi et al. (2003) found relatively little luminosity 
evolution in a study of 9000 early-type galaxies from the 
SDSS, in the redshift range 0 < z < 0.3. Furthermore, 
the metallicity evolution in star-forming galaxies from 
z ~ 0.1 to z ~ 0.4 is ~ 0.1 dex (Lara-Lopez et al. 2009). 
The metallicity evolution of early-type galaxies would 
not be expected to change more dramatically than this 
over our redshift range. Thus, considering the quoted 
uncertainty for Lboi.oiii of ±0.38 dex, we do not expect 
Eqn. 1 to be strongly affected by changing abundances 
in our redshift range up to z ~ 0.5. However, we will 
also check the significance of our results after excluding 
sources with z > 0.2. 

To estimate the jet mechanical luminosity from the 
1.4 GHz radio luminosity, we use the relation of Cav- 
agnolo et al. (2010), 

imech = 7.3 X 10 36 (L!.4 GHz/10 24 W Hz" 1 ) 0 ' 7 . (2) 

This scaling relation between the jet mechanical lumi¬ 
nosity (Lmech) and the 1.4 GHz radio luminosity is de¬ 
rived from X-ray cavity measurements in mainly low 
power radio galaxies with radio luminosities less than 
10 25 W Hz -1 (e.g. Birzan et al. 2008). The large scat¬ 
ter in the data (~0.8 dex) means that there is consider¬ 
able uncertainty in the slope of this relation, that makes 
its application to high-power radio sources questionable, 
as well as potentially introducing a systematic error in 
Tmech between the HERG and LERG samples because of 
their different luminosity distributions. However, God¬ 
frey & Shabala (2013) used an independent method for 
determining the total power of sources with radio lu¬ 
minosities greater than 10 25 W Hz -1 , specifically FR2- 
type sources. Despite the significantly different radiative- 
efficiencies of FR2 and FR1 type radio sources, Godfrey 
& Shabala (2013) surprisingly found very good agree¬ 
ment with the Cavagnolo et al. (2010) relation up to 
10 28 W Hz -1 . While this agreement is encouraging, the 
applicability of both the same normalisation and slope 
of Eqn. 2 across such a wide range in radio luminosity 
remains uncertain, with further work on a larger num¬ 
ber of sources across the full radio-jet luminosity range 
required. 

In order to calculate the Eddington luminosity, 

±Edd = 1-3 x 1O 31 M bh /M 0 W, (3) 

for each of our sources, we require an estimate of the 
black hole mass (Mbh)- Using the velocity dispersion 
of the host galaxy (cr*) taken from the SDSS spectrum, 
Mbh can be estimated for a large fraction of our sample 
through the M-cr* relation (e.g. Tremaine et al. 2002). 
This relation is based on the well determined dependence 
of the stellar velocity dispersion of a galaxy bulge on its 
black hole mass, described by 

log(M B H/M 0 ) = 8.13 ± 4.02 log(cr*/200 km s _1 ). (4) 

The SDSS stellar velocity dispersion estimates are ob¬ 
tained by fitting a spectral template across the rest-frame 
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wavelength range 400-700 nm, with emission lines regions 
explicitly masked out. The velocity dispersion is only es¬ 
timated for z < 0.4 and we have excluded any values 
above 420 km s -1 and below 70 km s _1 , since they can¬ 
not be reliably measured from the SDSS spectra. We also 
excluded any spectra with a median per-pixel signal-to- 
noise of less than 10. 14 This leads to 546 reliable veloc¬ 
ity dispersion values (22 HERGs and 524 LERGs) with a 
median redshift of 0.15 and with only four sources with 
z > 0.3. Thus, we find reliable velocity dispersion mea¬ 
surements for 40% of polarized HERGs and ~71% of po¬ 
larized LERGs. The lower fraction for HERGs is mainly 
due to a larger fraction of sources at higher redshift and 
with lower signal-to-noise. 

The stellar velocity dispersion as measured from the 
fixed 3” fiber aperture of the SDSS does not provide an 
accurate estimate of A7 bh in disk-dominated galaxies be¬ 
cause stars from the extended regions of the galaxy will 
be included. All the polarized HERGs and LERGs in 
our sample have a concentration-index (C = R90/R50) 
greater than 2.6, where the parameter i?90(i?50) rep¬ 
resents the radius that encloses 90% (50%) of the host 
galaxy optical light. A concentration-index greater than 
2.6 means that they are bulge dominated systems (Shi- 
masaku et al. 2001), as expected for the typical early- 
type hosts of radio-loud AGN. This means that the ve¬ 
locity dispersion as measured from the fixed 3” SDSS 
fibers (~8 kpc at z = 0.15) is representative of the bulge 
in the majority of cases. However, at redshifts, z > 0.2, 
the inclusion of the majority of the host galaxy light may 
affect the velocity dispersion measurement and the Mbh 
estimate in a systematic manner. At redshifts, 2 < 0.1 
and 0.1 < 2 < 0.2, we found median velocity dispersions 
of 240 km s _1 and 241 km s _1 , respectively, while from 
0.2 < z < 0.3, the median velocity dispersion increases 
to 257 km s -1 . This means that Eddington-scaled ac¬ 
cretion rates for sources with z > 0.2 are typically ~1.6 
times systematically larger than sources at z < 0.2. 

In Figure [l6j we plot Hi. 4 ghz versus the Eddington- 
scaled total accretion rate, corresponding to (L mec h + 
Tboi,oiii)/^Edd, for 22 polarised HERGs and 524 polar¬ 
ized LERGs. We see that the HERGs have systematically 
higher accretion rates (~ 0.5—10%) than the majority of 
the LERGs that have accretion rates <1%. In addition 
to this, it is clear that the more weakly accreting LERGs 
are more highly polarized in comparison to the HERGs 
and the strongly accreting LERGs. As noted by Mingo 
et al. (2014), the LERGs with Eddington ratios greater 
than ~1% may be mis-identified as radiatively inefficient 
sources and/or their radio luminosity is boosted by re¬ 
siding in dense environments. Thus, it is interesting to 
test the significance of the difference in Hi. 4 ghz between 
sources separated purely by their total Eddington ratios. 
A two-sided KS test indicates a highly significant differ¬ 
ence in Hi 4 GHz for sources separated at a total accretion 
rate of 0.5% of Eddington (p-value of 1.9 x 10 -15 , cor¬ 
responding to ~8cr for a normally-distributed process). 
We separated the sample at 0.5% of Eddington, instead 
of 1% of Eddington, in order to include all the HERG 
sources and better match the numbers in the two bins. 
If we exclude sources with z > 0.2, due to their system- 

14 See http://classic.sdss.org/dr2/algorithms/veldisp.html for 
details. 


atically higher Eddington luminosities, the KS test gives 
a p-value of 8.4 x 10 -9 (~6cr). 

Mingo et al. (2014) found that the mid-IR bolometric 
luminosity (Lboi jr) provided a more reliable estimate of 
the accretion rate in HERGs. Thus, we use this as a 
consistency check against the accretion rates estimated 
using the [OIII] emission line. From the WISE 22/i?n 
flux density, we can estimate LboijR from the scaling 
relation of Runnoe et al. (2012), where log(Lboi,m) = 
15.035 + 0.688 log(XL\). While we find that Lboi.iR can 
be up to an order of magnitude larger than Lboi,om 
for the HERGs, the qualitative difference when Lboi.iR 
is used to estimate the total accretion rate instead of 
Tboi.oui does not result in a significantly cleaner separa¬ 
tion in the Eddington ratios between polarized HERGs 
and LERGs. This also does not hugely effect the signif¬ 
icance of the difference in Hi,4 GHz for sources separated 
at an Eddington ratio of 0.5% (p-value of 1.3 x 10 -8 , 
~5.7cr). 

The bolometric luminosity correction in Eqn. 1 is for 
the [OIII] flux uncorrected for dust extinction. This cor¬ 
rection can be important because while the [OIII] line 
emission region lies outside the dusty torus of the AGN 
(i.e. in the Narrow Line Region), it may still suffer sig¬ 
nificant amounts of extinction due to interstellar dust in 
the host galaxy. To investigate this potential error in 
Tboi.ouii we corrected for the effect of dust extinction 
on the [OIII] flux using the Balmer decrement (i.e. the 
ratio of the fluxes of the Ha and H/3 narrow emission 
lines, see Lara-Lopez et al. (2009) for details). The me¬ 
dian size of the dust correction to the [OIII] luminosity 
for our sample is ^0.2. Kauffmann & Heckman (2009) 
used X-ray, optical and mid-infrared spectra of a com¬ 
plete flux-limited sample of SDSS type 2 AGN to derive 
a mean bolometric correction of 600 for the extinction- 
corrected [OIII] luminosity. Using this bolometric correc¬ 
tion to the extinction-corrected [OIII] luminosity, instead 
of Eqn. 1, decreases the significance in ni4GHz between 
low and high Eddington ratios to ~4.8cr. However, this 
decrease in significance is unsurprising since we lose a 
large fraction of our sources because we can only obtain 
the extinction correction for ~50% of sources due to the 
absence of Ha and/or H/3 measurements. 

While the different bolometric correction factors pro¬ 
duce slightly different results and their application at 
redshifts z > 0.2 is questionable, the main result of 
this section, that Hi. 4 ghz can achieve higher values for 
sources with lower total Eddington accretion rates, re¬ 
mains valid, for a separation at ^0.5% of Eddington. 

4. DISCUSSION 

Our main observational result is the discovery of a 
difference in the integrated degree of polarization at 
1.4 GHz (Hi 4 GHz) between radio-loud radiative-mode 
AGN (i.e. HERGs and radio-loud QSOs) and jet-mode 
AGN (i.e. LERGs), with the weakly accreting sources 
extending to higher values of Hi. 4 ghz- We now discuss 
the potential significance of these results in linking the 
large-scale host galaxy environment with the accretion 
state of the supermassive black hole and the production 
of powerful radio-loud jets. 

4.1. Environments of HERGs and LERGs 
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At low redshift (z < 0.4), HERGs are found mainly 
in low galaxy-density environments (i.e. groups), while 
LERGs are found in a wide range of environments from 
poor to rich groups and clusters (e.g. Best 2004, Hardcas- 
tle 2004). They also differ in their host galaxy types, with 
HERGs mainly residing in galaxies with significant on¬ 
going star formation implying an abundance of cold gas, 
with LERGs mainly hosted by massive elliptical galaxies 
with hot X-ray halos and a minimal supply of cold gas. 
This has led to the development of models claiming that 
the two different types of AGN result from their presence 
in different environments. 

In particular, Hardcastlc et al. (2007) proposed that 
the difference between HERGs and LERGs is due to the 
source of the accreting gas: Bondi accretion of the hot, 
X-ray gas in ellipticals for LERGs, and cold gas accretion 
from disk-galaxies required for HERGs. Alternatively, 
the type of radio-loud AGN activity may instead be trig¬ 
gered by a fundamental change in the accretion disk 
structure at an Eddington-scaled accretion rate of ~1% 
(e.g. Mingo et al. 2014), with a radiatively-inefficient ac¬ 
cretion flow below, and a radiatively-efficent thin disk 
above, ~1% of Eddington. However, these two scenarios 
are not mutually exclusive and both could potentially 
act together. It remains unclear as to what the exact 
nature of the relationship is between the large scale envi¬ 
ronmental properties of AGN and the type of accretion 
onto their supermassive black holes. 

In Section 3.5, we have shown that our sample of radio- 
polarized AGN display similar host galaxy properties and 
accretion rates as the full sample in BH12 and elsewhere 
(i.e. the host galaxies of our sample of HERGs typically 
have higher rates of star formation, lower black hole mass 
and higher accretion rates than the LERGs). Thus, we 
can use our sample of polarized sources to reliably in¬ 
fer some of the key properties of the HERG and LERG 
populations in general. 

The polarized sources in our sample are in less dense 
environments (in terms of the galaxy number density 
within 0.5 Mpc) than unpolarised sources. Furthermore, 
it is the most underdense environments of the polar¬ 
ized LERG sources that typically have higher values of 
ni.4GHz (Section 3.4.3). With a similar approach, Shi et 
al. (2010) studied the most highly polarized sources in 
the NVSS and found a median galaxy count within a 1 
Mpc radius of 7 ± 3 (for only eight sources) compared to 


9 ± 7 for a much larger low-polarization control sample. 
Their result is not inconsistent with our findings, and is 
suggestive of an influence of galaxy number density on 
Hi .4 GHz■ Overall, galaxy number density estimates seem 
to be crude but statistically useful probes of the environ¬ 
ment densities of radio sources in relation to whether 
they are in a poor or rich group, or cluster. 

In considering the radio source environment as a key 
variable in our analysis, we also need to recognise that 
bent or disturbed morphology radio sources are often 
found in dense environments (e.g. Owen & Rudnick 
1976). In fact, radio sources with bent morphologies are 
used as a way to find candidate galaxy clusters (e.g. Wing 
& Blanton 2011). In Section 3.3.2, we showed that bent 
radio sources have significantly lower values of Hi. 4 ghz 
than straight ones. We also found that Hi. 4 ghz for 
LERG sources was weakly correlated with the total inten¬ 
sity spectral index (Section 3.4.4), providing additional 
evidence that the LERGs with low values of Hi. 4 ghz are 
preferentially located in denser environments. This re¬ 
duction in Hi 4 GHz could be due to a combination of 
polarization angle cancellation, due to asymmetric jet 
structure, as well as the large amounts of external Fara¬ 
day rotation expected in dense environments, that would 
cause strong depolarization of the radio emission. 

These results are generally consistent with Best (2004) 
who found, for z < 0.1, a strong correlation between ra¬ 
dio luminosity and environment richness for LERGs, but 
not for HERGs. Similarly, Ineson et al. (2013) found, 
for radio sources around z ~ 0.4, a weak correlation be¬ 
tween radio luminosity and cluster richness for LERGs, 
but no correlation for HERGs. However, if the environ¬ 
ment of radio sources does indeed strongly affect the in¬ 
tegrated degree of polarization, the question still remains 
as to how exactly this is achieved. Thus, a more sensi¬ 
tive probe of the immediate environment of radio sources 
is desirable for probing the link between galaxy environ¬ 
ment, black hole accretion and powerful jet production 
in more detail. One of the most sensitive probes of the 
local environment of radio sources is through Faraday ro¬ 
tation and its associated depolarisation. In Sections 3.4.1 
& 3.4.2, we showed that the currently available Faraday 
rotation and depolarisation data for our sources is not 
sufficient to probe the local source environment in de¬ 
tail, however, it is clear that the majority of depolarisa¬ 
tion occurs local to the source (e.g. Banfield et al. 2014, 
Fames et al. 2014). 

From our current investigation, we expect that the 
LERGs with high values of Hi. 4 ghz are suffering less de¬ 
polarization in less dense environments than the LERGs 
with low values of Hi. 4 ghz- This is supported by our ob¬ 
servation that sources with straight jets are much more 
likely than bent jets to have Hi. 4 ghz > 15% (Section 
3.3.2), and the environment density estimates discussed 
above indicate that the more highly polarized sources are 
more likely to exist in less dense environments. From 
simulations, the amount of depolarization is expected 
to decrease with increasing lobe size (e.g. Hardcastle & 
Krause 2014). However, we find that the sources with 
Lli. 4 ghz > 15% have smaller linear sizes (Section 3.3.3), 
possibly indicating that these LERGs also have intrin¬ 
sically weaker jets. While the differences between FR1 
and FR2 radio morphologies are generally attributed to 
a combination of jet power and how they interact with 
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their ambient environments (i.e. the less powerful jets 
of FRls are gradually decelerated by entrainment while 
the more powerful FR2s are not), our results in Section 
3.3.1 suggest that the FR morphology is not the primary 
driver causing the observed difference in II! 4 ghz of the 
LERGs and HERGs. 

Our claim that the density of the ambient environ¬ 
ment is the key variable for the spread in II1.4GHZ for 
LERGs does not explain the fact that HERGs are limited 
to ni.4GHz < 15%, since they are expected to exist in less 
dense environments, on average, than LERGs. A possi¬ 
ble remedy to this could be related to the much greater 
radiative output of the HERG central engine. The high 
radiative output is likely photo-ionising large amounts 
of the host galaxy interstellar medium and increasing 
the amount of ‘Faraday-active’ material (i.e. the amount 
of magnetoionic material contributing to the total Fara¬ 
day rotation as seen by the jet). The larger amount 
of magnetoionic material would also mean larger vari¬ 
ations in Faraday rotation causing greater amounts of 
depolarization in HERGs, on average, than in LERGs. 
There is a well studied correlation between the AGN 
accretion disk luminosity and the amount of the host 
galaxy that is ionised (e.g. Netzer et al. 2006). In fact, 
in some cases there is evidence of saturation of the ion¬ 
isation radius, indicating that the AGN has potentially 
ionised the entire galaxy (e.g. Curran & Whiting 2010, 
Hainline et al. 2014). However, the central AGN need 
not be responsible for the ionisation of the large scale 
host galaxy environment in all cases. Extended emission 
line regions (EELRs), on scales of ~100 kpc, are often 
found in sources with strong nuclear emission lines and 
are coincident with regions of high depolarization in the 
lobes of radio galaxies (e.g. Pedelty et al. 1989). There 
is also statistical evidence that these EELRs are aligned 
with the radio source axis and have similar extents to 
the radio lobes (McCarthy et al. 1987, Baum & Heck¬ 
man 1989). This suggests that the radio jet may also 
be responsible for the ionisation of its immediate en¬ 
vironment through shocks, and shock-induced star for¬ 
mation, as well as entrainment. Interestingly, from a 
Faraday rotation measure study of 26 radio-loud AGN 
(0.3 < z < 1.3), Goodlet & Kaiser (2005) found that the 
difference in the RM of each lobe and the dispersion in 
RM both correlate with redshift, suggesting more dense 
magnetoionic environments at higher redshift. Since all 
their high-redshift sources are HERGs and their low- 
redshift sources are dominated by LERGs, this lends sup¬ 
port to the argument that HERGs maximum integrated 
polarization may be limited due to large amounts of Fara¬ 
day depolarization local to the source. Once higher pre¬ 
cision Faraday rotation measure and depolarisation data 
are obtained for a large number of HERGs and LERGs, 
it would then be interesting to consider other factors, in 
addition to the magnetoionic environment, such as the 
frequency of gas-rich mergers as a function of redshift 
and environment, the star formation histories of the host 
galaxies and their X-ray properties, for example. 

Banfield et al. (2014) also studied the 1.4 GHz po¬ 
larization properties of a large sample of radio galaxies 
and radio-loud QSOs taken from the catalog presented 
in Hammond et al. (2012). Banfield et al. (2014) dis¬ 
cussed the difference in the degree of polarization be¬ 
tween radio galaxies and radio-loud QSOs as a result 


of a cosmic evolution in the space density of quiescent 
galaxies (i.e. the typical hosts of LERGs). They also 
found no significant difference in the RM distribution 
of the two types of sources but concluded that most of 
the extragalactic RM must originate close to the source. 
Such studies have important implications for predict¬ 
ing the number density of polarized sources that will 
be detected in deep polarization surveys with the SKA 
and its precursors/pathfinders. The integrated degree 
of polarization of radio-loud AGN at 1.4 GHz has been 
claimed to be anti-correlated with the total radio flux 
density (e.g. Mesa et al. 2002, Tucci et al. 2004, Taylor 
et al. 2007, Subrahmanyan et al. 2010, Stil et al. 2014), 
and although this relation is disputed (Hales et al. 2014), 
it was claimed that it is the intrinsically less luminous ra¬ 
dio sources that are more highly polarized (Banfield et 
al. 2011, 2014). Preliminary efforts were made to explain 
the anti-correlation based on a change in the dominant 
population of polarized radio sources from FR2s to FRls 
towards lower flux densities (O’Sullivan et al. 2008, Stil 
et al. 2014). However, the results of this paper strongly 
suggest this is not the case. We also find a weak but sig¬ 
nificant anti-correlation between the 1.4 GHz radio lumi¬ 
nosity and Hi 4 GHz (Section 3.2) and instead claim that 
it is most likely driven by the difference in Hi. 4 ghz be¬ 
tween radio-loud radiative-mode AGN, that dominate at 
high luminosities, and the jet-mode AGN, that dominate 
at lower luminosities. 

4.2. Jet production efficiency 

After liking the large scale magnetoionic environment 
of radio sources to their host galaxy accretion states 
(i.e. HERGs and LERGs), it is interesting to discuss the 
potential implications of the accretion of magnetised gas 
from large distances in enabling the production of pow¬ 
erful, radio-loud jets. 

The existence of large amounts of poloidal magnetic 
flux close to a supermassive black hole is synonymous 
with jet production (e.g. McKinney & Blandford 2009, 
Porth et al. 2011, Fendt et al. 2014) and also with the 
extraction of the rotational energy of the black hole in 
the so called magnetically-arrested accretion disk regime 
(MAD; Igumenshchev et al. 2003). In MAD systems, it 
is expected that the mass accretion rate and jet magnetic 
flux is strongly correlated (Tchekovskoy et al. 2011), 
which is supported by recent observational results (Za- 
maninasab et al. 2014, Zdziarski et al. 2014). However, 
these observations are of sources with radiatively-efficient 
accretion and the MAD systems have only been success¬ 
fully simulated in radiatively-inefficient regimes. Fur¬ 
thermore, van Velzen et al. (2013) have argued that the 
MAD model is inconsistent with the radiatively-efficient 
accretion in radio-loud QSOs, on the grounds that the 
spread in the optical-radio luminosity relation does not 
allow for a large range of black hole spins. Mocz & Guo 
(2015) presented a model for interpreting MADs in both 
the radiatively-efficient and radiatively-inefficient accre¬ 
tion regimes, however this is based on the assumption 
a MAD can be successfully formed in the radiatively- 
efficient accretion regime. 

In any case, the required amounts of magnetic flux to 
generate such strong poloidal fields close to the black hole 
most likely needs to be advected from large distances, 
since the generation of such large magnetic flux in-situ is 
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considered unlikely (e.g. Begelman 2014). If the differ¬ 
ences in accretion mode of HERGs and LERGs is mainly 
due to the large scale environment (i.e. HERGs have a 
radiatively-efficient accretion flow due to a large supply 
of cold gas and LERGs have a radiatively-inefficient ac¬ 
cretion flow due to a hot interstellar medium), then the 
ability of the accretion flow to form a MAD system and 
produce powerful radio-loud jets may be related to the 
large scale magnetised environment. Thus, the study of 
the large scale magnetoionic environment of HERGs and 
LERGs through their polarization and Faraday rotation 
properties may provide important constraints on the like¬ 
lihood of formation of a MAD. 

In Section 3.5.3, we found that the sources with low 
Eddington-scaled accretion rates (Ltotai/TEdd < 0.5%) 
could achieve higher values of Hi. 4 gh z than the sources 
with Ltotai/iEdd > 0.5%. This could be explained 
by a more uniform jet/lobe magnetic field structure 
and/or less Faraday depolarization in the environments 
of sources with Ltotai/AEdd < 0.5%. Our analysis of the 
effect of radio morphology on Hi. 4 ghz suggests that in¬ 
trinsic magnetic held differences are not the dominant 
effect. We also find no significant difference between the 
RMs of strongly and weakly accreting sources, but we 
claim that the currently available RM data are dom¬ 
inated by the RM contribution from the Milky Way 
and thus insufficient to determine the true RM contri¬ 
bution local to the source. Higher precision RMs from 
broadband and high angular resolution radio observa¬ 
tions should provide a robust test of this scenario through 
direct estimation of the amount of Faraday depolariza¬ 
tion local to the source. 

AGN that accrete in a radiatively efficient manner are 
considered to be the dominant type of AGN at high 
redshift and possibly responsible for the co-evolution of 
supermassive black holes and their host galaxies. Why 
only a small fraction (<10%) of these radiatively-efficient 
AGN are radio-loud remains an open question. One pro¬ 
posed solution is the ‘magnetic flux paradigm’ of Sikora 
& Begelman (2013), who propose that the radiatively- 
efficient radio-loud AGN must have undergone a previous 
radiatively-inefficient accretion phase in order to build up 
the large amounts of magnetic flux required for powerful 
jet production (c.f. Lubow et al. 1994). One way in which 
this could be achieved is by a major merger between a 
disk galaxy and a giant elliptical galaxy that has been un¬ 
dergoing hot accretion for some time. Such a scenario has 
some observational support in that the most luminous 
AGN are often triggered by mergers (e.g. Ramos Almeida 
et al. 2012). Whether or not giant elliptical galaxies have 
sufficient amounts of coherent magnetic fields that can 
be dragged inwards to the black hole horizon is currently 
unclear; however, this is not inconsistent with our obser¬ 
vations of highly polarized LERGs whose elliptical galaxy 
hosts would cause smaller amounts of Faraday depolar¬ 
ization if large-scale coherent magnetic fields were indeed 
present. Upcoming broadband radio polarization data 
from surveys, such as POSSUM on the ASKAP telescope 
(Gaensler et al. 2010) and the proposed VLA Sky Sur¬ 
vey, can provide important constraints on this proposed 
scenario. 

5. CONCLUSIONS 


We have presented a radio polarization study of 796 
radio-loud AGN (z < 0.7) in relation to their host galaxy 
properties, accretion states and large scale environments. 
The integrated degree of polarization at 1.4 GHz was ob¬ 
tained from the NVSS in the region that overlapped with 
the SDSS survey area. We used the SDSS optical spec¬ 
troscopic classifications of Best & Heckman (2012), who 
defined a large sample of NVSS radio sources into high- 
excitation and low-excitation radio galaxies (i.e. HERGs 
and LERGs). We find a fundamental difference in the po¬ 
larization properties between HERGs and LERGs, where 
the LERGs are observed to span the full range of ex¬ 
pected integrated degrees of polarization at 1.4 GHz 
(IIi.4 ghz), with a maximum of ~ 30%, while HERGs 
(and radio-loud QSOs) are restricted to a maximum of 
IIi.4 GHz ~ 15%. We also find that the weakly accret¬ 
ing LERG sources (with Eddington ratios < 0.5%) can 
attain higher values of Hi,4 GHz than the strongly accret¬ 
ing LERGs and HERGs (with Eddington ratios > 0.5%). 
This may be important for attempts to explain the pro¬ 
duction of powerful radio-loud jets by the ability of AGN 
host galaxies to accumulate large amounts of magnetic 
flux close to the black hole (e.g. Sikora & Begelman 
2013). Infrared data from WISE allowed us to deter¬ 
mine the host galaxy type, with LERGs mostly inhab¬ 
iting early-type elliptical galaxies and HERGs typically 
with disk-dominated hosts. This is broadly consistent 
with detailed optical studies which find that the host 
galaxies of HERGs are typically of lower stellar mass and 
have bluer colours than LERG hosts (e.g. Heckman & 
Best 2014). Thus, we maintain that the results from our 
sample of polarized HERGs and LERGs are applicable 
to the HERG and LERG populations in general. 

The scenario that we consider to best explain the cur¬ 
rent data is one in which the full range of polarization 
from < 1% to ~30% of the LERGs, can be explained 
by the full range of environments in which LERGs are 
observed to reside in, in the local universe. In high gas 
density environments, the jet emission structure is more 
likely to be bent, causing polarization angle cancellation, 
while significant Faraday depolarization caused by the 
high density magnetoionic medium can help to further 
suppress the observed integrated polarization at 1.4 GHz. 
In poor group environments, the jets are more likely 
to remain straight and along with the lower expected 
amount of external Faraday depolarization, this leads 
to the highest integrated degree of polarization. Addi¬ 
tionally, the LERG sources with n 14 GHz > 15% have 
significantly smaller linear sizes than the other LERGs, 
potentially indicating that these LERGs may also have 
intrinsically weaker jets. The restricted integrated de¬ 
gree of polarization in HERGs is not easily explained 
since the HERGs reside in lower galaxy density environ¬ 
ment, on average, than the LERGs. However, we suggest 
that the high ionising luminosity of the central engine 
ionises a significant fraction of the host galaxy, increas¬ 
ing the amount of magnetoionic material contributing to 
the Faraday rotation local to the source. This would be 
expected to generate large amounts of Faraday depolar¬ 
ization in the HERG environment, potentially restricting 
the integrated degree of polarization at 1.4 GHz from 
reaching values greater than ~15%. 

However, although we expect Faraday effects to domi¬ 
nate at 1.4 GHz, we do not find any direct evidence for 
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this from the Faraday rotation measures (RMs) derived 
from the NVSS data in Taylor et al. (2009). These RMs 
have been shown to be dominated by the magnetoionic 
material of the Milky Way and are relatively insensitive 
probes of the local environments of extragalactic radio 
sources (e.g. Oppermann et al. 2015), thus, we cannot 
test our above hypothesis in a robust manner. We have 
also considered intrinsic magnetic field differences of the 
large scale radio morphological classes (i.e. FR1, FR2 and 
FRO) as a potential origin of the difference in the polar¬ 
ization properties of the HERGs and LERGs. However, 
we do not find any conclusive evidence suggesting that 
the FR morphological class is the primary driver. To 
understand the relative importance of Faraday rotation 
and intrinsic magnetic field structure in determining the 
polarization properties of HERGs and LERGs requires 
high precision RM and polarization measurement across 
a wide range of frequencies. Upcoming broadband radio 
polarization surveys, such as the Polarization Sky Sur¬ 
vey of the Universes Magnetism (POSSUM) on the Aus¬ 
tralian Square Kilometre Array (SKA) Pathfinder tele¬ 
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APPENDIX 

THE BRIGHTEST RADIO SOURCES: I > 100 mjy 

Stil et al. (2014) showed that the median degree of polarization of all NVSS sources increases towards lower flux 
densities (e.g. the presented fit to their data gives a median degree of polarization of 1.93% at 100 mjy and 2.17% at 10 
mjy). Therefore, we split our sample at I = 100 mjy in order to compare the properties of the brighter sources to the 
fainter sources in our sample. At I > 100 mjy, there are 211 LERGs and 37 HERGs. We find no significant difference 
in III.4 ghz between HERGs and LERGs for I > 100 mjy. The number of LERGs in our sample with Hi. 4 ghz > 15% 
decreases from 71 at I < 100 mjy (~10% of LERGs) to only 2 at I > 100 mjy (~1% of LERGs). 

In Figure Al, we repeat Figures 4, 6, 7 & 8 from the main text but now only include sources with I > 100 mjy. 
Comparing Fig. Al(a) to Fig. 4, we see that there is no strong redshift dependence on flux in our sample. The median 
redshift for LERGs changes from 0.18 to 0.16, and from 0.25 to 0.23 for HERGs. The majority (88%) of the radio-loud 
QSOs have I > 100 mjy. In Fig. Al(b) & (c), we see that most sources with Li. 4 ghz < 10 24 W Hz -1 are removed. In 
fact, the median radio luminosity increases by a factor of approximately 2 for both HERGs and LERGs and there is 
no longer any significant evidence for a correlation between Hi. 4 ghz and Li. 4 GHz- Fig. Al(d) confirms that even for 
the brightest sources, the radio morphology does not show any clear split in radio luminosity into FR1 and FR2 type 
sources. There is no significant difference in Hi. 4 ghz between FR1, FR2 or FRO source morphologies. 

Interestingly, the median total intensity spectral index of sources with I > 100 mjy is —0.7 while it is —0.6 for 
I < 100 mjy. If the spectral index is a useful probe of the environment density (as argued in Section 3.4.4), then the 
less-steep spectral index for sources with / < 100 mjy may explain why there is a greater fraction of LERGs with 
Lli.4 ghz > 15% at these lower fluxes. This provides some support for our conclusion that a sub-population of polarized 
LERGs with weaker jets and smaller sizes in under dense environments can achieve Hi. 4 ghz > 15%. We have spectral 
index measurements for ~60% of sources with / > 100 mjy, but for only ~20% of sources with I < 100 mjy. The 
upcoming broadband radio polarization surveys will provide much better statistics in order to investigate this in more 
detail. 
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Fig. Al.— Figures (a), (b), (c) & (d) correspond to Figures 4, 6, 7 & 8 of the main text, respectively, but now only including sources 
with total radio flux density greater than 100 mjy. 





























